Anisotropic Strong Coupling Calculation of the Local Electromagnetic Response of 

High-Tc Superconductors 

A. Bille and K. Scharnberg 
Fachbereich Physik, Universitdt Hamburg, Jungiusstrafle 11, D-20355 Hamburg, Germany 

The electromagnetic response of the Cu02-planes is calculated within a strong coupling theory 
using model tight binding bands and momentum dependent pairing interactions representing spin 
fluctuation and phonon exchange. The superconducting state resulting from these interactions has 
d-wave symmetry. With phonon exchange included the order parameter amplitude grows rapidly 
below Tc at elevated frequencies which leads to improved agreement with the observed temperature 
dependence of the penetration depth. Good agreement between theory and experiment can only be 
achieved if it is assumed that the strength of the quasiparticle interaction decreases with temperature 
in the superconducting state. The amount of this reduction depends sensitively on the momentum 
dependence of the interactions, the energy dispersion and the position of the Fermi line. 



The fundamental information required to perform 
strong coupling calculations for conventional supercon- 
ductors is the Eliashberg function which incorporates 
both the phonon spectrum and the clcctron-phonon inter- 
action. These properties can be determined in the normal 
state and they are not expected to to be affected signifi- 
cantly by the onset of superconductivity. Since much ev- 
idence for a highly anisotropic pair state in the high-Tc 
materials has accumulated in the recent past, the mo- 
mentum dependence of the Eliashberg function needs to 
be taken into account. Upon the assumption that the un- 
derlying interaction is an exchange of spin- fluctuations, 
such an Eliashberg function has been modelled quite suc- 
cessfully to explain normal state properties. There is still 
some controversy about the best possible model but we 
feel that the function suggested by Monthoux and Pinesjd 
which contains the momentum dependence in a nonsep- 
arable form, is a reasonable starting point. Since this 
Eliashberg function is supposed to arise from electron- 
electron interactions, we expect it to change when super- 
conductivity sets in. One way this change could manifest 
itself is a lower frequency cut-off approximately of the 
size and temperature dependence of the superconducting 
order parameter amplitude.H 

In this paper we have not included such strongly tem- 
perature dependent modifications of the Eliashberg func- 
tion. Instead, we focus on the momentum and frequency 
dependence of the various self-energies, in particular the 
scattering rate, which in conjunction with the anisotropy 
of the Fermi velocity could have a marked effect an the 
conductivity. For an energy dispersion of the formu 

e(k) — —t [2 cos{kxTr) + 2 cos(/cj,7r) 

—AB cos{kxTT) cos{kyTr) + /i] , (1) 

with kx,y given in units of tt/o, we have plotted in Fig. 
1 the Fermi lines (inset) and the Fermi velocities in the 
first quadrant of the Brillouin zone. For each choice of 
the parameter B the band filling parameter /i is cho- 
sen such that the Fermi line meets the zone boundary 
at the same point (1,0.1). For a simple tight binding 



band B — 0, quasiparticles with momenta in (1, 1) di- 
rection have indeed the highest Fermi velocities and thus 
contribute most strongly to the conductivity. For the 
modified tight binding band with B = 0.45 there is very 
little variation of the Fermi velocity, which may be an 
artefact of this particular model band structure, as the 
comparison with the dot-dot-dashed curves show. 
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FIG. 1. Fermi velocities |Ve(k)p = (hn/a)'^ as calcu- 
lated from Eq. ^ for a band width 8t = 1.44 eV. The inset 
shows the corresponding Fermi lines. The dot-dot-dashed 
lines are for the 8 parameter tight binding approximiation to 
the odd band of the Cu02-double layers in YBCO,u scaled 
to the same band width. 

It has been noted before that a pairing interaction, 
which is attractive at small momentum transfers, can con- 
tribute to the formation of a d-wave stateJj In addition 
to the spin fiuctuation induced pairing interactioru 
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we have, therefore incorporated into our strong coupling 
calculations a phonon contribution of the form 



a'^Fphiq, n) = gph -G(q) V 

TT < * 
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Numerical results are for = 41 meV and F = 8 meV. 
A very similar frequency dependence for a'^Fp^ lias been 
suggested by GolubovE Note that at low frequencies both 
functions vary linear Iv-, with frequency but for a typical 
value ujsF = 7.7meVp which is much smaller than the 
value 30meV used in Ref. [ ^, o^Fsf dominates. 

We have antisymmetrized a'^Fph{^.^^) with respect to 
f2, which has little effect on the physical meaning of this 
function but ensures that the coupling constant 
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Xph^2 / dq^ / dqy / -j^-a'^ Fph{q,fl) = gph (4) 
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is independent of both Hq and F.B The normalized q- 
dependent form factor 



G(q) = 36.9 6-25?' 
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reflects nesting properties of the Fermi surface. The ex- 
pression Eq. (^) represents a good fit to the numerical 
results shown in Fig. 3 of Ref. ||. In all our calcula- 
tions the coupling constant 175^ was adjusted such that 
Tc = 92 K was kept constant. 

Details of the calculation of the self-energies Z(k,a;), 
which renormalizes the frequency to a) = LuZCkjUj), 
x(k, w), which is added to the band energy Eq. (p, and 
of the anomalous self-energy <i>(k,a;)la will have to be 
given elsewhere. To calculate these quantities we used 
the Eliashberg equations in the form given by Marsiglio 
et al. [Ref. |lof . 

For the energy dispersion Eq. with B = 0.45 we 
show in Fig. 2 ReZ(k,a;), which is related to the mass 
renormalization, as function of uj for three different tem- 
peratures and two points on the Fermi line, one on the 
zone boundary, the other one at the position of the order 
parameter node. In Fig. 3 the corresponding results are 
shown for Inict;Z(k, w), which, is often identified with the 
quasiparticle scattering ratdiil. 

These figures show that there is considerable 
anisotropy already in the normal state, indicating much 
stronger interactions where the Fermi line meets the zone 
boundary than on the zone diagonal. This is due to the 
fact that for the particular band considered the nesting 
condition q = (±1,±1) which maximizes the interaction 
Eq. (H), is fulfilled to a much higher degree for points 
kF = (1.0,0.12) than for kj. = (0.37,0.37). Large val- 
ues of Z are found on the image of the Fermi line undes 
a translation by Q, i.e. near k = (1 — 0.37, 1 — 0.37).B 
The anisotropy of Z on the Fermi line is not altered when 
phonons are included because these scatter most strongly 
between Fermi lines near ki? = (1.0, ±0.12). An interac- 
tion peaked at zero momentum transfer, however, would 
reduce the ansitropy. 

Both mass renormalization and scattering rates are 
smaller at low frequencies when phonons are included 



because the coupling parameter Xsfi defined in anal- 
ogy to Xph Eq. (|), has to be taken as Xsf — 2.545 eV 
to obtain Tc — 92 K which is larger than Asf + Xph — 
1.979 eV + 0.25 eV which gives the same Tc- The source 
of this discrepancy is the contribution to the phonon spec- 
trum at high frequencies which is absent in o^Fsf- 

When the temperature is lowered below Tc, the scat- 
tering rate drops much faster at the zone boundary where 
it was largest. The effect that this behaviour will have 
on the conductivity depends very much on the variation 
of the Fermi velocity (see Fig. 1). At T « O.lTc, the 
scattering rate for uj < 15meV has drojj>P|ed well below 
typical disorder induced scattering rates.EJ. 
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FIG. 2. Mass renormalization function versus frequency 
for reduced temperatures T/Tc = 1.0, 0.9, 0.1. At 
Lj — 60 meV, Re Z increases monotonically with decreasing 
temperature. Dashed curves are for Qph = in Eq. , solid 
curves are for gph ~ 0.25 eV. Note the different ordinate 
scales. 
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FIG. 3. Scattering rates versus frequency for 
T/Tc = 1.0,0.9,0.1. At around to = 30meV, ujlmZ in- 
creases monotonically with increasing temperature. Dashed 
curves: Qph = 0, solid curves: Qph ~ 0.25 eV. 

In Fig. 4 we show the order parameter A(k, w) = 
<i>(k, w)/Z(k, w), on the Fermi line at the zone bound- 
ary. The ratio Re A(w = 0,T = 0.1Tc) /Tc = 2.79 without 
phonons and 3.42 with phonons is larger than the weak 
coupling value 2.14 but is still at the lower end of the 
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range of values rmjiired to fit experiments witliin weak- 
coupling theories. El 
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FIG. 4. Real and imaginary parts of the order parame- 
ter versus frequency for T/Tc — 0.99, 0.9, O.f and gph = 
(dashed) and Qph = 0.25 eV (solid). 

The surface impedance is calculated in the local limit 
which is well justified for all high-Tc materials. The con- 
ductivity is obtained without taking vertex corrections 
into account. Since the interactiops are momentum de- 
pendent, this is an approximation.113 




FIG. 5. Real part of the conductivity versus reduced tem- 
perature at fixed microwave frequency. 

Fig. 5 shows the real part ai of the conductivity at 
87 GHz as function of temperature for different energy 
dispersions and different interactions. The only parame- 
ter adjusted in this figure is gsF and this is chosen to fit 
Tc, not (Ti. For B — 0.45 in Eq. (|l|) the agreement with 
a typical experiment (Ref. [ |Tl|, Sample B) is good, be- 
cause the scattering rates near the nodes are much smaller 
than the average value. At low temperatures, ai drops 
too fast because the inelastic scattering rate falls below 
the experimental frequency. For these temperatures it is 
essential to take disorder induced elastic scattering into 
account. This would also reduce the peak height so that 
the conclusion seems to be inescapable that the interac- 



tions must become smaller below Tc- How large a reduc- 
tion one would infer from experiment depends crucially 
on the energy dispersion. With B = we have nearly 
perfect nesting and hence a large scattering rate on the 
Fermi line where the order parameter vanishes and where 
the Fermi velocity is maximal. For this case we find that 
the peak in cti is far too small to be anywhere near the 
data. 
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FIG. 6. Imaginary part of the conductivity 

Fig. 6 compares the superfluid density in the form 
a2{uJ,T) with data taken at 87GHz.Eil Since ReZ renor- 
malizes the bare London penetration depth, which is de- 
termined by the band structure, we obtain A(0) = 123 nm 
when phonon-exchange is included instead of A(0) = 
155 nm. The data are scaled to fit these different val- 
ues of A(0). Including phonons gives a much improved fit 
because of the fast rise in ReA(ki?,ti;) at high frequen- 
cies (see Fig. 4). As in the case of cti, reaching agreement 
with experiment seems to require scattering rates which 
drop more rapidly in the superconducting state than our 
temperature-independent Eliashberg fmictions predict. 
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